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ABSTRACT 

Equations for the wave-averaged three-dimensional momentum equations have been published in 
this journal. It appears that these equations are not consistent with the known depth-integrated 
momentum balance, especially over a sloping bottom. These equations should thus be considered 
with caution as they can produce erroneous flows, in particular outside of the surf zone. It is 
suggested that the inconsistency in the equations may arise from the different averaging operators 
applied to the different terms of the momentum equation. It is concluded that other forms of the 
momentum equations, expressed in terms of the quasi-Eulcrian velocity, are better suited for three 
dimensional modelling of wave-current interactions. 



1. Introduction 

The wave-averaged conservation of momentum can take 
essentially two forms, one for the mean flow momentum 
only, and the alternative form for the full momentum, which 
includes the wav e pseudo-mome ntum (hereinafter 'wave 
momentum', sec lMcIntvrei[l98lh . This question is well 
known for depth-integrated equat ions ([Longuet-Higgins and 
Il964t lGarrettl[T976l : ISmith|[2006h . but the vertical profiles 
of the mass and moment um bal a nces are more complex. 
The pioneering effort of iMellorl (|2003l hereinafter M03) 
produced practical wave-averaged for the total momentum 
that, in principle, may be used in primitive equation mod- 
els to investigate coast al flows, suc h as t he wave-driven 
circulatio ns observed b y iLentz et aL ( 2008 ). The first for- 
mulation (|Mellorll2003h was slightly inconsistent due to the 
improper approximation of wave motion with Airy wave 
theory, which is not enough on a sloping bottom, however 
small the slop e may be. This question was discussed by 
lArdhuin et all (|2008t) . and a correction was given and ver- 
ified. These authors acknowledged that these equations, 
when using the proper approximation, are not well suited 
for practical applications because very complex wave mod- 
els are required for the correct estimation of the vertical 
fluxes of wave momentum, that are part of the fluxes of 
total momentum. 

Although M03 gave correct wave-forcing expressions - 
in terms of velocity, pressure an d wave-induc ed displace- 
ment, before any approximation - IMellorl (I2008L hereinafter 
M08) derived a new and different solution from scratch. 



The two theories may be consistent over a flat bottom, but 
they differ at their lowest order over sloping bottoms, so 
that the M08 equ ations are likely to be flawed, given the 
analysis of M03 bv lArdhuin et al. ( 20081 ). and the fact that 
their consistency was not verified numerically over sloping 
bottoms. 

instead, M08 asserted that the equat ions are consis- 



— — — msteaa, iviuo asserted tnat tne equat ions are consis 
^eitf with the depth- integrated equations of Phillips dl977t) 



Furth er, about the test case proposed by lArdhuin et al 



(|2008h . M08 stated that the wave energy was unchanged 
along the wave propagation and that the resulting wave 
forcing whould be uniform over the depth. Here we show 
that the M08 equati ons do not y ield the known depth- 
integrated equations (jPhillipsllf 9771 ) with a difference that 
produces very different mean sea level variations when waves 
propagate o yer a sloping b o ttom. As for the test case of 
proposed bv lArdhuin et al. ( 2008 ). we show that a consis- 
tent analysis should take into account the small but sig- 
nificant change in wave energy due to shoaling. In the 
absence of dissipative processes, the M08 equations can 
produce spurious velocities of at least 30 cm/s, with I m 
high waves over a bottom slope of the order of 1% in 4 m 
water depth. 

2. Depth-integration of the M08 equations 

For simplicity we consider motions limited to a verti- 
cal plane (x, z) with constant water density and no Cori- 
olis force nor wind stress or bottom friction. The wave- 
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averaged momentum equation in M08 takes the form 

dU dU 2 dUW 
dt dx dz 
and the continuity equation is 

dU dW 
dx dz 



= -Sf + F 
ox 



0. 



(1) 



(2) 



Where U and W are the Lagrangian mean velocity com- 
ponents, which contains the current and Stokes drift veloc- 
ities, g is the acceleration due to gravity and rj is the time- 
averaged water level at the horizontal position x. The force 
given by M08 on the right hand side of ((T|) can be written 
as the sum 

F = F"°*+F nn 



of a wave-induced pressure gradient 



F, 



Mi is 



dS. 
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and the divergences of the horizontal flux of wave momen- 
tum, 



F, 



:Uus 



dS uu 



du 2 
dx 



dx 

— (kEFocFcs) 



(6) 
(7) 



where E is the wave energy, k is the wavenumber, u and w 
are respectively the horizontal and vertical wave-induced 
(orbital) velocities. Ejj is defined by 



Ed = if z ^ rj and 



Endz = — . 



(8) 



Using the mean water depth D, and bottom elevation 
—h, Feci Fss and F$c are non-dimensional functions of 
kz and kD, 



Fee 
Fss 



F 



sc — 



cosh(fcz + kh) 
cosh(kD) 

sinh(fcz + kh) 
s'mh(kD) ' 

sinh(fcz + kh) 



cosh(kD) 

The depth-averaged mass-transport velocity is 
1 rv 



U 

M08 correctly noted that 



D 



Udz. 
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(9) 
(10) 

(11) 
(12) 
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with S^x 7 given bv lPhillipsI (|l977l ). However, for a depth- 
uniform t7, th e depth integrated momentum equation in 
Phillips! (Il977h is 



-gD 



drj d 



: P77 



(14) 



The forcing in the depth-integration of (JTJ) differs from 
the forcing in (fT4|). because the gradient is inside of the 
integral, namely, 



rj q cA/08 

* X dz 
u dx 



dx 
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(Z = -h) 



dh 
dx 



(15) 



The depth integral of M08 thus includes two extra term. In 
particular S%, m (z = -h)dh/dx = -2kE(dh/dx)/ sinh(2fc£)) 
can be dominant over a sloping bottom. As a result the 
momentum balance in M08, unlike M03, does not produce 
the known set-down and set-up. This is illustra t ed in figure 
[1] We take the case proposed bv lArdhuin et al.l (|2008l ) with 
steady monochromatic waves shoaling on a slope without 
breaking nor bottom friction and for an inviscid fluid, con- 
ditions in which exact numerical solutions are known. The 
bottom slopes smoothly from a depth D = 6 to D = 4 m. 
We added a symmetric slope back down to 6 m to allow 
periodic boundary conditions if needed. For a wave pe- 
riod of 5.24 s the group velocity varies little from 4.89 to 
4.64 m s _1 , giving a 2.7% increase of wave amplitude on 
the shoal. Contrary to statements in M08, dE/dx is sig- 
nificant, with a 5.4% change of E over a few wavel engths. 

F rom the Eulerian analysis of that situation (e.g. lLonguet-Higgins 
19671 ). the mean water level should be 0.32 mm lower on 



the shoal (figure [T]). iRivero and Arcilla (Il995h established 
that there is no other dynamical effect: a steady Eulerian 
mean current develops, compensating for th e divergence of 
the w ave- induced mass transport (see also lArdhuin et al 
20081 ). 



3. Flows produced by the M08 equations 

Because the relative variation in phase speed is impor- 
tant, from 6.54 to 5.65 m s _1 , the Stokes drift accelerates 
on the shoal. The Eulerian velocity u is irrotational, thus 
nearly depth-uniform, and compensates the Stokes drift 
divergence by a convergence. The Lagrangian velocity U, 
shown in figure [2l is the sum of the two steady velocity 
fields. 

We now solve for the equations derived by M08. The 
numerical sol ution is obtaine d by coupling the WAVEWATCH III 
wave model (jTolmanl [2009h . solving the phase-averaged 
wave action equation, and the MAR S3D ocean circula- 
tion model ( Lazure and Dumas] 2008). This c oupling uses 
the generic coupler PALM (jBuis et al.l l2008h . The feed- 
back from flow to waves is negligible here and was thus 
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Fig. 1. Mean sea surface elevation induced by monochromatic waves propagating over the smooth bottom shown in 
figure [21 with amplitude H s = 0.34 m and period T = 5.24 s. The extra terms forcing terms in eq. (fl"5|) lead to an 
ov erestimation of the se t-down by more than 50 % for this case. ARB08 stands for quasi-Eulerian momentum equations 
of lArdhuin et all (|2008t ). 



turned off. MARS 3D was implemented with 100 sigma 
levels regularly spaced, and 5 active points in the transver- 
sal y direction, with 2 extra wall points, and 2 ghost points 
needed to define finite differences, it is thus a reall three- 
dimensional calculation although the physical situation is 
two-dimensional. There are 78 active points in the x di- 
rection. The time step was set to 0.05 s for tests with 
H s = 1.02 m (1 s for H s = 0.34 m). For simplicity, the 
wave model forcing is updated at each time step. We use 
Eq. ([T]) transformed to ? coordinates, with <; defined by 
z = s(x,q,t) = rj + cP + s (|Melloill200a ). 



dU 
dt 



U 



r dU 
dx 



WdU 87} 
zT7k~ ~ 9 dx~' 



(16) 



where the advection terms are obtained by using Eq. ([2]). 

The flow boundary conditions are open. The monochro- 
matic wave amplitude a = 0.12 m translates into a signifi- 
cant wave height H s of 0.34 m for random waves with the 



same energy. We also test the model with a = 0.36 m, i.e. 
H s = 1.02 m, still far from the breaking limit in 4 m depth. 

The discontinuity of the vertical profile in the forc- 
ing F, due to the Ed term, is not easily ingested by the 
numerical model, and generates a strongly oscillating ve- 
locity profile (Fig. [3]). These oscillations are absent at 
depths larger than 0.8 m, consistent with the zero val- 
ues of F below the surface. A realistic constant viscosity 
K z — 2.8. 10~ 3 m 2 s -1 removes the oscillations by diffus- 
ing the negative term —dEp/dx over the vertical. Yet this 
term is a momentum source that produces velocities one 
order of magnitude larger than the Stokes drift U s , with 
an opposite sign (Fig. [3J . The spurious velocities given by 
M08 with a realistic mixing are most pronounced for waves 
in not too shallow water (Table [1} , and comparable with 
those given by the M03 equations without mixing. 
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Fig. 2. Lagrangian velocity U for the inviscid sloping bottom case with H s = 1.02 m and T = 5.24 s, obtained from the 
quasi-Eulerian analysis as U = u + U s . Contours are equally spaced from -0.01 to 0.025 m s _1 . The thick black line is 
the bottom elevation. 



Table 1. Model results with Mellor (2008): Sur- 
face velocity at x =200 m (on hte up-slopc) for different 
model settings . The settings corresponding to the test in 
Ardhuin et al.l ( 2008 ) are given in the second line. The 
surface velocity values are written for the time t = 900 s 
except for the case without mixing (t = 360 s). 



H s (m) 


T p (b) 


K z ( 


m 2 a 1 ) 


U (m s- 1 ) 


1.02 


5.6 







0.6116 


0.34 


5.6 







0.2127 


0.34 


13 







0.3164 


1.02 


5.6 


2.8 


.10~ 3 


-0.1594 


0.34 


5.6 


2.8 


.10~ 3 


-0.0256 


0.34 


13 


2.8 


TO" 3 


-0.0007 



4. Conclusions 

We showed that the equations derived bv lMellorl (|2008h 
appear inconsistent with the know depth-integrated mo- 
mentum balances in the presence of a sloping bottom. In 
the absence of dissipation, a numerical integration of these 
equation produce unrealistic surface elevations and cur- 
rents. The currents may reach significant values for very 
moderate waves, exceeding the expected results by one or- 
der of magnitude. While we d id not discuss the origin of the 
inconsistency, it appears that lMellor (12008ft used a different 
averaging for the pressure gradient term and for the advec- 



tion terms of the same equation. We believe that this is the 
original reason for the problems discussed here. The spu- 
rious velocities produced by M08 are likely to be dwarfed 
by the strong forcing imposed by breaking waves in the 
surf zone. Nevertheless, we expect that the M08 equations 
can produce large errors for continental shelf applications, 
such as the investigation of cross-shore transports outside 
of the surf zone. Alternativel y, equations for th e quas i- 



Eulerian veloc i ty ca n be used [M cWil 



jams et al.l (|2004f ) 



Ardhuin et all (|2008[ ); IUchivama et "all (|2009l) which do not 
have such problems. 
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